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Measurement of the photon structure function at ALEPH 

A. B6hrer a Si-99-03 

a Fachbereich Physik, University of Siegen, 
57068 Siegen, Germany 

The photon structure function F% has been measured with data taken by the ALEPH collaboration at LEP 
centre-of-mass energies y/s = 91 GeV with (Q 2 ) of 9.9, 20.7 and 284 GeV 2 and y/s = 183 GeV with (Q 2 ) of 13.7 
and 56.5 GeV 2 . For the data at y/s = 183 GeV a two-dimensional unfolding method employing the principle of 
maximum entropy is used, which reduces the errors compared to one- dimensional methods. 



1. INTRODUCTION 

The photon structure function in e + e~ col- 
lisions is measured [0J^] where one of the incident 
beam leptons is scattered at sufficiently large an- 
gle to be detected. In these singly-tagged events 
the negative momentum transfer Q 2 = —q 2 of 
the photon emitted by the tagged lepton is Q 2 = 
2i?ta g £ ; beam(l — cos 9) . The second lepton stays 
undetected so the momentum transfer is small. 

The process then can be viewed as inelas- 
tic electron-photon scattering, where a quasi-real 
photon is probed by a virtual photon B : 
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The inelasticity, measured with y, is small and 
the term with F^ may be neglected. Equation [l] 
can be used to relate the distribution in x and Q 2 
to the structure function F% . A more elaborate 
introduction and discussion of F% you may find 
inRef. g. 

In this article the recent measurements of F2 
with the ALEPH experiment are presented 
Measurements at two different centre-of-mass en- 
ergies and various Q 2 have been performed. The 



data are used to test different model parameteri- 
zations of F^ . Special emphasis has been put on 
a two-dimensional unfolding method to extract 
the true x distribution dN/ dx from the measured 
data. 



2. DATA 

The ALEPH detector and its performance have 
been described in detail elsewhere |J. Charged 
tracks and neutral calorimeter energy as defined 
by the ALEPH energy flow package |^| are used 
in these analyses. 

Single tag events are selected by the lepton 
detected in the electromagnetic calorimeters of 
ALEPH: ECAL and LCAL for yfi = 91 GeV 
(LCAL and SiCAL for y/s = 183 GeV). The scat- 
tering angle is measured in the range of —0.6 < 
cos6> tag < 0.95 in ECAL and 65mrad < tag < 
150mrad in LCAL (60mrad < 8 tag < 155mrad 
at 183 GeV). Further cuts on E tag and a veto on 
a second tag are applied. 

The visible hadronic final state has to con- 
sist of at least three charged tracks and has to 
have an invariant mass W v - ls > 2 GeV. Fur- 
ther cuts for rejection of beam-gas events are ap- 
plied. Additional cuts are required for the data 
at =91 GeV rejecting background from Z de- 
cays. 

We are left with clean data samples with a re- 
maining background of a few percent. The sam- 
ples are divided into three (two) subsamples; see 
Table |. 

The data are corrected for trigger efficiency 
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Table 1 

Data samples with their statistics and Q 2 





#of 
events 


Q 2 range 
GeV 2 


(Q ) 
GeV 2 


E 

^cms 

GeV 


ECAL 


163 


35 - 3000 


284 


91 


LCAL 


1647 


13 - 44 


20.67 


91 


LCAL 


1543 


6 - 13 


9.93 


91 


SiCAL 
LCAL 


1208 
861 


7- 24 
17 - 200 


13.7 

56.5 


183 
183 



(which is close to 100%), downscaling, and ac- 
ceptance. For the correction it is important to 
note that the invariant mass W and therefore x 
is poorly measured only, see Fig. 0. This is es- 
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Figure 1. x v - ls versus Xtmc 
events at t/s = 91 GeV 



for LCAL tagged 



pecially due to the significant portion of the en- 
ergy of the event which goes in the forward re- 
gion where there is little or no tracking informa- 
tion. Therefore extraction of a measurement of 
the photon structure function from the data de- 
pends on the use of a model of the production of 
hadronic final states from 7*7 collisions. There 
is at present no complete theoretical description 



of this process, so a number of phenomenological 
models are used. 

Models giving reasonable descriptions of global 
event variables are then used to correct the data 
using an unfolding method. 

The models chosen at y/s = 91 GeV are firstly a 
combination of QPM and VDM HQ , both with 
JETSET fragmentation ij, and secondly HER- 
WIG § with GRV LO |L0). QPM and VDM 
are combined to form a single set, weighting each 
sample so as to give the best x 2 between distri- 
butions (the number of energy flow objects, the 
transverse momentum of all energy flow objects 
with respect to the beam direction, and the thrust 
of the event) predicted by the combined simu- 
lation and the data. In the HERWIG program 
the defaults for version 5.9 are used, apart from 
the intrinsic transverse momentum fct of the par- 
tons in the target photon. This was modified ac- 
cording to the scheme proposed by ZEUS and de- 
scribed in Ref. 

The models chosen at y/s = 183 GeV are firstly 
HERWIG § with GRV LO @, secondly HER- 
WIG § with SaS |]| and thirdly PHOJET @ 
(only used for SiCAL data, because the descrip- 
tion of the LCAL data is only moderate). 



3. UNFOLDING AND 
EXTRACTION OF F 2 7 

As pointed out in the previous section, the in- 
variant mass W is poorly measured by current 
detectors and an unfolding method has to be em- 
ployed. Both at y/s = 91 GeV and at 183 GeV a 
regularization procedure is used to suppress os- 
cillations in the result. 

At = 91 GeV the unfolding was performed 
using the procedure proposed by Blobel JjlJ. 
This procedure fits a sum of spline curves to the 
data after passing them through the x v - la versus 
£true response matrix obtained from the simu- 
lated events (HERWIG and QPM+VDM), sup- 
pressing oscillations which have higher frequency. 
The structure function Fj is then obtained using 
the G ALUG A program |ll| , where is set to its 
theoretical value. With F^ set to 1 the output is 
used as reference for the extraction of F^ of the 
data. 
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At y/s = 183 GeV a new, recently proposed 
method jll] has been used for the unfolding 
(Ref. H and references therein). It uses the 
principle of maximum entropy. In addition, this 
method allows a two-dimensional unfolding to be 
applied. The regularization function S(£t) used is 
the Shannon entropy fllTf: 



M , , 

^(M) = -E— ) 

Mtot VMtot/ 



(4) 



when are the parameters to be estimated in 
bin j and /Ztot is their sum. This entropy-based 
regularization makes no reference to the relative 
locations of any of the bins. Therefore the princi- 
ple of maximum entropy can be directly applied 
to multidimensional distributions. This method 
reduces the model dependence by including for 
each event not only x, but in addition some other 
variable characterizing the final state. For this, 
the variable En has been introduced, defined as 
the total energy of the particles having angles 
with respect to the beam line less than 17°. Here 
the x resolution degrades considerably for increas- 
ing En, see Fig. [| 
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Figure 2. a; V i S versus Xtruc f° r subsequent bins of 
the variable En using HERWIG simulation. 



The two-dimensional unfolding leads to smaller 
statistical errors. This is because in the one- 
dimensional case, the effective weight of each 
event is determined by the average x resolution in 
the bin. With two-dimensional unfolding, those 
events with low En are given a higher weight 
in the final result. The improvements achievable 
by two-dimensional unfolding where investigated 
quantitatively in Ref. Jilt] - ^ n example is shown 
in Fig. H where a toy-MC sample was used and 
two different response matrices to demonstrate 
the differences in one- and two-dimensional un- 
folding. Nevertheless the model dependence re- 
mains a significant source of systematic uncer- 
tainty, and several event generators are used in 
the presented analysis to account for this. After 
unfolding the structure function F% is obtained 
from the x distribution dN/dx using the same 
MC simulation that was used for the determina- 
tion of the response matrix. The measurement 
in data is taken as the average of the distribu- 
tions obtained with the response matrices of the 
models used. 
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Figure 3. One- (upper) and two-dimensional 
(lower plots) test unfolding in SiCAL- (left) and 
LCAL- (right) tagged events using a sample of 
toy-MC. 
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4. RESULTS 

4.1. Measurements for F% 

The final results for F^f are shown in Fig. |] for 
V^i = 91 GeV and in Fig. | for ^ = 183 GeV, 
resp. The inner error bars are statistical errors, 
the total error bars represent the quadratic sum 
of statistical and systematic errors. The system- 
atic error comprise uncertainties from the unfold- 
ing (spread of the results based on the different 
models etc.), from detector efficiency and accep- 
tance including trigger, selection criteria etc. The 
91 GeV analysis includes a small theoretical error 
from the assumed form factor of the virtual pho- 
ton as well as the dependence on the assumption 
oiFl. 
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Figure 5. The values for F^/a at y/s = 183 GeV 
((Q 2 ) =13.7 and 56.5 GeV 2 ) compared to three 
different parameterizations. 



4.2. Comparison of F^ to parameteriza- 
tions 

For both centre-of-mass energies and the differ- 
ent Q 2 ranges all of the parameterizations provide 
a good description of the data for x > 0.1. For 
lower x the LAC values are significantly too high. 
At y/s — 91 GeV a test has been made using event 
reweighting to ensure that the difference is not an 
artifact of the unfolding procedure; for more de- 
tails see Ref. Q. The low x region is sensitive 
to the gluon content of the photon. A x 2 was 
calculated in the eleven bins in x to quantify the 
comparison to various sets of the parton density 
function of the photon, see Table 0. Those that 
show significantly large values of x such as LAC 
1 and 2 and Whit 4,5, and 6, contain a large gluon 
content, resulting in a rapid rise in the structure 
function at low x. 



Figure 4. The values for F^ja. at yfs = 91 GeV 
compared to three different parameterizations. 



4.3. Dependence of F^ from (Q 2 ) 

The values of F^ for the three Q 2 ranges at 
y/s = 91 GeV have been averaged over the region 
0.1 < x < 0.6. A logarithmic rise with (Q 2 ) is 
seen as expected from theoretical predictions us- 
ing the parameterizations listed in Table 0. 
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Table 2 

The values of x 2 obtained data at to the more re- 
cently calculated photon parton density function, 
see [0 and references therein. 
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5. CONCLUSIONS 

Single-tagged two-photon events recorded by 
the ALEPH detector at LEP I and LEP II have 
been studied in three and two bins of Q 2 , which 
have a mean of 9.9, 20.7 and 284 GeV 2 and 13.7 
and 56.5 GeV 2 . The data have been used to mea- 
sure the hadronic structure function F% as a func- 
tion of x. The comparison with parameterizations 
show that those with parton density functions 
that contain a large gluon content are inconsis- 
tent with data. The rise of with (Q 2 ) is found 
compatible with the available parameterizations. 

A two-dimensional unfolding technique using 
the principle of maximum entropy has been suc- 
cessfully applied; with the second variable chosen 
as En, defined as the total energy of the parti- 
cles having angles with respect to the beam line 
less than 17°. This unfolding method reduces 
the statistical errors and the model dependence of 
the extracted as compared to one-dimensional 
procedures. 
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